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INVESTIGATION OF AN OPTIMUM DETECTION SCHEME 
FOR A STAR-FIELD MAPPING SYSTEM* 

By Melvin D. Aldridge and Leonard Credeur 
Langley Research Center 

SUMMARY 

An investigation was made to determine the optimum detection scheme for a star- 
field mapping system that uses coded detection resulting from starlight shining through 
specially arranged multiple slits of a reticle. With threshold levels established by the 
Neyman-Pearson criterion, the computer solutions of equations derived from a theoreti- 
cal model are used to judge performance, and then practical problems such as synchro- 
nization, code construction, and detector noise are considered. 

On the basis of the Neyman-Pearson criterion it was determined that the maximum 
probability of detection occurred with the use of a single transparent slit in the modulating 
reticle. However, use of multiple slits improved the system's ability to reject the detec- 
tion of undesirable lower intensity stars, but only by decreasing the detection probability 
of the lower intensity stars to be mapped. Also, it was found that the coding arrangement 
affected the root-mean-square star-position error and that detection is possible with 
error in the system's detected spin rate, though at a reduced probability. 

INTRODUCTION 

The problem investigated herein was encountered at the Langley Research Center 
during the design of a star -field mapping system for a spinning rocket probe (ref. 1). The 
system uses a coded detection scheme resulting from starlight shining through specially 
arranged multiple slits of a reticle. It appears desirable to use a number of slits so that 
more than one pulse per star will be available for a cross-correlation detection process 
with the known slit pattern. 

The purpose of this study is to determine the optimum number and arrangement of 
slits (i.e., optimum codes) for maximum probability of detecting stars in the star scanner. 

*Some of the material presented herein formed part of a thesis entitled "Investiga- 
tion of the Optimum Detection System for a Star Field Mapping System," offered by 
Melvin D. Aldridge in partial fulfillment of the requirements for the degree of Master of 
Electrical Engineering, University of Virginia, Charlottesville, Virginia, August 1965. 



In addition, the ability of the system to reject the detection of undesirable lower intensity 
stars is considered. Performance information is obtained by utilizing computer solutions 
of equations derived from a theoretical model and then considering practical problems 
such as synchronization, code construction, and detector noise. 

Although star mapping provided the prime motivation for this study, the results can 
in general be applied to a light-source detection system of the type considered. 

SYMBOLS 

aj decision that xj is true state of affairs 

1 + B secondary emission noise factor 

b additional noise factor 

Cj cost of falsely detecting a star 

Cjj cost of deciding that j occurred when i really occurred 

C m cost of missing a star 

D equivalent dark current rate, photoelectrons per sample 

e electron charge, coulombs 

f r code repetition bit rate 

G average photomultiplier gain 

Ij 3 average photomultiplier anode dark current, amperes 

i,j events 

k average number of photoelectrons emitted per second 

k n average noise photoelectron rate per reticle slit 

1 number of times 1 is detected as 1 
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m number of times 0 is detected as 1 

N number of transparent slits in modulating reticle 

n number of photoelectrons 

P D probability of detecting a star, assuming code is fully in correlator 

P D l’ P D 2 P D ^ or cross_corr el a tion detection methods 1 and 2, respectively 
P D,6 probability of detection for code in correlation position 6 

PpA probability of false alarm 

P FA,1> P FA,2 P FA ^ or cross_co:r relation detection methods 1 and 2, respectively 
P FA,G given (allowable) Pp^ 

P(g,N) probability that g events will occur out of N possible events 

P(i) probability that event i will occur 

P(i/j) probability that event i will occur, knowing that event j has occurred 

Pjj probability that i is detected as j (i,j either 0 or 1) 

P(n) probability that n photoelectrons will be emitted during time r 

Pp probability that star is present for detection 

q number of times 0 is detected as 0 

R overall average risk of decision scheme 

Tq average cost of error when star is not present 

r j average cost of error when star is present 

s integral multiple of exact code repetition rate 
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Tg bit decision threshold 

Tc correlation decision threshold 

r ^'C,l>^’C,2 f° r cross-correlation detection methods 1 and 2, respectively 

t time 

At uncertainty in sample time 

tg star pulse width in time (bit time) 

t 0 beginning of detection code 

t s sample length in time, tg/s 

condition i is true state of affairs (i either 0 or 1) 
y result of an observation 

a variable of integration of standard Gaussian distribution 

r signal-to-noise ratio per reticle slit (ratio of average number of star photo- 

electrons to average noise photoelectrons emitted per second) 

y average gain per dynode 

6 offset of detected code in correlator 

e rms position error of star 

X average noise photoelectrons per sample per slit 

v error factor of sampling repetition rate 

a T 2 total variance of theoretical statistics 

r integration or counting time, seconds 

u> vehicle spin rate, radians/second 

4 



Subscripts: 


max maximum 

min minimum 


GENERAL DESCRIPTION OF SYSTEM 

The star -field mapping system for a spin-stabilized rocket probe which is con- 
sidered in this report can be represented by figure 1. Light from the star to be detected 
is concentrated on the cathode of a photomultiplier tube by a telescope that has a field of 
view directed out of the side of the probe and is normal to the spin axis. As the vehicle 
rotates about its spin axis, this field of view sweeps a band in the celestial sphere. An 
opaque reticle with transparent slits parallel to the spin axis is positioned in the focal 
plane of the optics. The spin of the probe causes star images to sweep across the trans- 
parent slits, resulting in intermittent exposure of the photocathode to the images. 

Thus, a train of coded current pulses appears at the tube's anode with average 
amplitude determined by the star's magnitude and format determined by the arrangement 
of the transparent slits. The form of these pulses after ideal detection is illustrated in 
figure 2. This coded pulse train will consist of two basic states: a 1 state occurs if a 
star is shining through any one of the slits, while a 0 state occurs if a star is not shining 
through a slit. If the coded pulse train is interpreted as a binary code, the l's and 0's 
are binary digits which are commonly called bits. 



Figure 1.- Basic star-mapping system. 
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Figure 2.- Ideal ‘bit-detection output signal. 

The 0 state before detection will have an average level determined by the background 
light from other stars and celestial bodies in the field of view. Thus, the 1 state before 
detection will have an average amplitude determined not only by the star magnitude, but 
also by the aforementioned background light. This background light is considered as 
noise since it hinders the detection of the star. Also contained in the pulse train before 
detection will be noise generated within the photomultiplier tube and photon noise due to 
the quantum nature of light. 

In reality, background light consists of discrete stars distributed over a wide range 
of magnitudes; however, it is convenient to consider that the aggregate of the very dim 
stars and celestial bodies constitutes a spatially homogeneous background light. Conse- 
quently, stars of sufficient brightness to allow reasonable detection probability can be 
considered as point light sources on top of the background. In the practical system it is 
desirable to map only stars whose intensity is sufficient to allow reliable detection, while 
rejecting (not detecting) all stars of less intensity. This will be referred to as the rejec- 
tion property of a detection scheme. 

SIMPLIFIED THEORETICAL ANALYSIS 

Because of the many dependent variables involved, it is expedient to conduct the 
theoretical analysis on a simplified basis. This section first presents such a simplifica- 
tion and the associated assumptions. After development of the equations necessary to 
describe the entire detection process, the results of computer solutions are presented. 

Simplification and Assumptions 

The system can be represented by the simplified block diagram in figure 3. Since 
the photomultiplier tube is a photon- or quantum-sensitive device, the basic detection 
process can be simplified by counting the discrete photoelectrons (see appendix for justi- 
fication). As shown in the block diagram, two separate decision processes take place in 
the overall detection system. First is the bit detection, which attempts to determine if a 
star is or is not shining through a slit. Second is the cross -correlation detection, which 
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determines whether an adequate number (with the proper arrangement or coding) of these 
bits have been detected to ascertain the presence of a star. We can suppose that the cor- 
relator consists of a shift register into which the results of the bit detection are serially 
fed. The contents of the register are compared at each bit time with the known code 
determined by the reticle slit configuration. This comparison is the basis of the cross- 
correlation detection process. 


Background 

light 



Star present Star not present 


Figure 5*- System block diagram. 


The following assumptions will be made for the simplified theoretical analysis: 

1. The bit detector counts discrete photoelectrons. 

2. Samples are taken only when the star image lies wholly within or not within a 
slit (i.e., bit synchronization is known). 

3. Either a complete code due to the presence of a star or only detected noise bits 
exist in the correlator. That is, the correlation that takes place while serially shifting 
the code bits (due to the presence of a star) into and out of the correlator will be neglected 
for now. 

4. Noise generated within the photomultiplier will be ignored. 

5. The effect of having two stars of sufficient intensity for reliable detection within 
the field of view simultaneously will be neglected. 

6. The average number of emitted photoelectrons per sampling period will be large 
enough to allow the Poisson distribution to be approximated by the normal distribution 
(ref. 2). 

7. The number of l's and the number of 0's for a given coding arrangement are 
taken to be equal. Also, the pulse widths (bit times) of the l's and 0's are assumed to be 
equal. 


7 


In addition, it is assumed that the resolution of the optics is much better than the 
width of a transparent slit. Thus the resulting pulse will have nearly instantaneous rise 
and decay times. Otherwise the problem involves time-variable statistics, which is 
beyond the scope of this paper. The effects on a practical system when assumptions 2, 3, 
and 4 are not made will be considered later in the paper. 

Bit-Detection Analysis 

Bit detection involves the counting of discrete photoelectrons which are emitted 
from the photocathode during the sample period. These emitted electrons will possess 
conventional shot noise and thus will have a Poisson statistical distribution (ref. 3). This 
distribution is expressed by the relation 

P(n) = ^ 7 - exp(-kr) ( 1 ) 

n.' 


where 

n number of photoelectrons 

k average number of photoelectrons emitted per second 

t counting or integration time, seconds 

P(n) probability that n photoelectrons will be emitted during time r 

The statistics can be completely described from a knowledge of the average signal 
and noise photoelectron rates and the integration time r. These rates are defined by 

k n average noise photoelectron rate per reticle slit 

r signal-to-noise ratio per reticle slit (ratio of average number of star photo- 

electrons to average noise photoelectrons emitted per second) 

If the background light is considered as having homogeneous spatial density over the 
entire celestial sphere, the total background noise falling on the photocathode becomes a 
linear function of the number of slits in the reticle. Since synchronization information 
is known by assumption 2 , r can be made equal to tg, where tg is the star pulse 
width in time (bit time). The total average number of photoelectrons emitted per sample 
period is 
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where 


Signal + Noise = k n t B (N + r) = X(N + r) 


( 2 ) 


N 


total number of transparent slits 


average noise photoelectrons per sample per slit, k n t B 


These relations can now be applied to the Poisson distribution to describe the bit 
detection statistics. A bit threshold T B establishes the decision level. If the count 
during an integration period falls below this level, the decision is made that a 0 was 
present, or that no star was shining through a slit. Similarly, if the count appears at T B 
or above, the decision is made that a 1 was present, or a star was shining through one of 
the slits. The probability that a 0 will be detected as a 0 is 

P 00= ^ ^-exp(-NX) (3) 

n=0 

and the probability that a 1 will be detected as a 1 is 


'11 


V [\(N + T)1 r -1 

= 2 1 exp[-X(N + r)J 

n=T R 


(4) 


As discussed by Feller (ref. 2), the Poisson distribution can be approximated with 
good accuracy by a Gaussian distribution whenever the average value NX is sufficiently 
large, as stated in assumption 6. Thus, in the rest of this paper equations (3) and (4) will 
be approximated by 


t b -nx 


Pqo = Yi \ 1 2NA exp(-a2)da 


(5) 


P ll=-F 



T B -x(N+D 

j2x(N+r) 


exp(-Q!2)da! 


( 6 ) 


As will be found later, these Gaussian forms are convenient for the inclusion of photomul- 
tiplier noise in the analysis. These relations express the probabilities of correct bit 
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decision. The probabilities of false bit decision are then 


P 01 = 1 " P 00 


( 7 ) 


p 10 = i - p ll 


( 8 ) 


where 

P 01 probability a 0 is detected as a 1 

P 1q probability a 1 is detected as a 0 

It should be noted that for the special case of only one slit, the bit detection and cor- 
relation detection are synonymous, and equations (3) to (8) completely describe the 
system. For this case the probability of detecting a star when it is present is Pjj and 
the probability of a false alarm is Pqj. 

Cross-Correlation Detection Analysis 

As a result of the bit-detection process, the coded signal for each star consists of a 
series of l's and 0’s whose signal format is referred to as PCM (pulse-code-modulated) 
nonreturn-to-zero format. Cross -correlation between the detected bit code and a known 
code, determined by the geometric pattern of the transparent slits in the reticle, forms 
the basis for the cross-correlation detection. Knowing the detection criterion, it is pos- 
sible to calculate the probability of detecting a star if it is assumed that all the code bits 
from the presence of a star have been shifted into the proper correlation position (see 
assumption 3). Similarly, the probability of false alarm can be calculated by assuming 
that the detected bit code in the correlator should consist of all 0 bits in the absence of a 
star. For analysis, two different cross-correlation schemes will be utilized. Many 
schemes are possible, allowing different weights to the proper detection of 0's. The two 
methods chosen are extremes of a sort: (1) No weight placed on proper 0 detection; 

(2) equal weight placed on proper 0 and 1 detection. 

Method 1 — "one " detection only.- The first technique is concerned only with the 1 
bits and does not consider the 0's that occur as a result of the opaque area between slits 
of the modulating reticle. A correlation threshold T^ j can be defined so that T^ j 
or more 1 bits must agree with the 1 bits of the known code for the presence of the star 
to be detected. The binominal statistical distribution applies, since there exist N pos- 
sible events, each with the same probability of occurrence. This distribution is described 
by the relation (ref. 4) 
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P(g,N)= - P(i) e [l ^ P(i)] N ' 6 (9) 

g.(N - g): L J 

where 

N total number of possible times event i can occur (number of slits) 

P(i) probability that event i will occur 

P(g,N) probability that event i will occur g times out of N possible times 


Thus the overall probabilities of star detection and false alarm for method 1 are 
expressed as 


l = T C,l 


N -l 


( 10 ) 


P FA,1 


I p o”( 1 - p oi) N ' m 

m=T C,l 


( 11 ) 


Method 2 — "zero -one" detection .- The second technique treats the proper detection 
of 0's and l's as being equally important. At this point it should be recalled that in 
assumption 7 both the pulse width and the number of the 0's are equal to those of the l’s. 

A threshold Tq 2 can defined so that the total sum of all properly detected l's and 0's 
must be equal to or greater than T^ ^ f° r the presence of a star to be determined. The 
value of T C .2 shall be greater than N. Although the statistics again follow from equa- 
tion (9), all possible combinations of 0's and l's that will give Tc,2 or greater must be 
included. Assuming that the 0's and l’s are independent, the probability of star detection 
for method 2 becomes 


N 


P D;2 - p(rs + O’B i t c _ 2 ) = 2 

;=t C;2 -ni 


N 


£ p 00< q ’ N) p ii(*> N ) 

q=T C,2 _ * 


( 12 ) 


where 


p Oota’ N) ■ P 0 ° q ( 1 - P Oo) 


N-q 
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Plltt,N) 


n: 

z:(n - 1 ): 



Similarly, for the probability of false alarm, 

N 

£ P 00 (q,N) P 0 l(m,N) 

q=T c>2 - m 

where 


FA, 2 


N 

I 

m=T C,2" N 


P 0 i(m, N ) 


N* p rn A 

m.'(N - m)J 01 V 



N-m 


( 13 ) 


Choice of Thresholds 

Although bit and correlation thresholds were defined in the preceding sections, 
their magnitudes were not discussed. Their proper choice depends directly on the estab- 
lishment of an optimum system criterion. One such criterion exists in Bayes’ decision 
rule, which defines the conditions that minimize the overall cost of detection errors 
(ref. 5). Threshold magnitudes will now be considered by applying Bayes’ decision rule. 

Let the true state of affairs be represented by xg when a star is not present and 
by xj when a star is present. Let y represent the results of an observation related 
to the true state of affairs; it can be related to the true state of affairs by two probability 
density functions in the following matrix form: 


y 

x o p (y/ x o) 

x i p (y/ x i) 


(14) 


If aj represents the decision that xj is the true state of affairs when xj is really the 
true state, the relationship can be represented by 


a 0 a l 


p ( a o/ x o) 

P ( a l/ x o) 

P ( a o/ x l) 

p( a l/ x l) 


A relative cost C^j 


of the above decision can be assigned as 


(15a) 
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(15b) 


. 

a 0 

a i 


a 0 

a l 

x 0 

n 1 

o 

o 

C 01 

x 0 

0 

Cf" 

X 1 

_ c io 

Cll 

= X 1 

c m 

0 


where 

c 00 = Cii =0 since no cost is associated with a true decision 
C 01 = C f cost of falsely deciding a star was present 

Cio = c m cost of missing a star 

The average cost of error when the star is present is 

r l = (16) 

The average cost of error when the star is actually not present is 

r o = CfP(ai/ Xo ) (17) 

Note that 

P^xjJ = Pp (probability that a star is present) 

p ( x o) = 1 - Pp 

P^ai/x^ = Pjj (probability of detecting a star) 

p ( a o/ X l) = 1 ‘ P D 

P( a l/x(^ = PpA (probability of false alarm) 

p ( a o/ x o) = 1 " P FA 

The overall average risk is 

R = Ppr j + (l - Pj)r 0 = P p C m (l - P^ + (l - Pp)c f P FA (18) 

From Bayes' decision rule, the threshold should be chosen so that R is minimized. 
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At this point it would be helpful to discuss the relative magnitudes of Pp, C m , 
and Cf for a practical mapping system. A typical field of view might be 6° with slits of 
approximately 1 minute of arc. This is only a small percentage of the total steradians 
scanned during a vehicle rotation. Because of the relatively low density of stars bright 
enough for reliable detection, only a small percentage of the total time of a vehicle's 
rotation is involved in a positive detection process. In other words, the system will be 
vulnerable to false alarms much longer than it will be involved in detecting a star. 
Therefore 

(l “ Pp) » Pp (19) 

If it is assumed that the cost of one false alarm is about the same as the cost of missing 
a star, Cf and Cm are of the same order of magnitude. Furthermore, it is always 
desirable to have 

P D > P FA ( 2 °) 

With the above relative magnitudes in mind, it can be seen from equation (18) that 
if the ratio Ppy/PpA I s max i m i ze( t, the total risk R is minimized. However, the max- 
imization of this ratio relies only on the relative rates of change of Pp and for 

a varying threshold. Since decreasing Pp also decreases Pp^.? anc * vice versa > such 
a maximization can occur for false rates much lower (or higher) than those allowable, 
which could give a Pp much lower than would have been obtained if the allowable PpA 
had been utilized. Thus, it appears convenient to establish an allowable PpA an( * then 
choose thresholds to minimize 1 - Pp or maximize Pp. This is known as the Neyman- 
Pearson criterion. 


Method of Computer Analysis 

Equations (5), (6), and (10) to (13) were programed on a digital computer. The 
Neyman-Pearson criterion was used by choosing an allowable false alarm probability 
PFA,G- A correlation threshold Tc was chosen, and Poo was determined from equa- 
tion (11) or (13) by an iterative process which gave PpA within 1 percent of PpA,G on 
the low side. From equation (5) the Tg which gave this Pqo to within 0.00001 percent 
on the low side was determined. The probabilities Pf j and Pp were then calculated. 
This procedure was followed for both methods of cross -correlation. By varying the cor- 
relation threshold Tq for given values of X, T, and N, an optimum correlation 
threshold could be found that gave the maximum possible Pp. Then by varying N and 
using its corresponding optimum Tq values for given values of X, r, and PpA G’ the 
variation of Pp was studied. This procedure was repeated for various combinations of 
X, r, and PpA,G. 
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Discussion of Theoretical System Analysis 

As indicated in the preceding section, the problem in calculating the maximum Pj) 
for given values of X , r, and N is the proper choice of a correlation threshold Tq. 
In figure 4, J?jy is shown as a function of Tq for various values of r. These curves 
show that an optimum threshold exists. However, note that as r increases the curve 
flattens, and thus the criticalness of the optimum choice decreases. Although it is not 
easily shown in graphical form, as T increases the optimum threshold level tends to 
decrease. These characteristics are true also for a change of any variable which 
increases Pjj (i.e., increasing Pjj increases Pj) and flattens the curve) . 
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(a) Cross-correlation method. 1: detection of l*s only. 



(b) Cross-correlation method 2: detection of l*s and 0 r s. 

Figure 4-.- Probability of detection as a function of correlation 
threshold. N = 7$ A = 100; %a,G = 
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The optimum number of slits for a maximization of the probability of detection will 
be considered next. In each of the following cases, for a given number of slits N, the 
corresponding optimum Tq was used. 

Figure 5 shows that increasing X increases P D . At first thought this might 
appear incorrect, since an increase in X involves increases in both background noise 
photoelectrons and star photoelectrons, and these are increased in the same proportion 
because r is constant. Inspection of the first and second moments of the Poisson dis- 
tribution shows that the mean increases as X, while the rms deviation increases as \fx. 



(a) Cross-^correlation method 1: detection of 

l T s only. 



(h) Cross-correlation method 2: detection 

of 1 1 s and 0 's . 


Figure 5-- Probability of detection as a function of 
number of slits and noise strength. F = 0-7; 

P FA,G = 10-5 * 



16 




Thus the ability to discriminate between the noise distribution, with mean AN, and the 
signal-plus-rioise distribution, with mean A(N + r), has improved although the noise 
level has increased. 

As expected, figures 6 and 7 show that increasing r or Pfa increases Pp. 

As N increases, PD decreases faster for lower r, and also Pd, 2 appears to 
decrease slightly faster than Pd,1 f° r given values of A, r, and PpA,G‘ These 
properties are better exhibited by plotting Pd as a function of T, as shown in figure 8 
for N = 1, 3, and 8. Since decreasing star brightness is analogous to decreasing r, 



(a) Cross-correlation method 1: detection 

of l's only. 


1.0 

.8 

.6 

P D, 2 

.4 

.2 

0 



(b) Cross-correlation method 2: detection 

of 1 * s and 0 1 s . 


Figure 6.- Probability of detection as a function of 
number of slits and signal-to-noise ratio. 

A = 100; P PA>G = 10-5. 
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these characteristics are the desirable ones for improving the rejection property, that is, 
the ability of the system to reject (not detect) stars of undesirable lesser intensity. Thus, 
increasing N increases the rejection property. For values of r of about 0.8 or 
higher, Pj) in figure 8 is almost the same for all values of N, but for lower r values, 
PD is considerably less for greater N. 

It is interesting that the optimum correlation threshold T q discussed earlier is 
used to the advantage of the code's rejection property. This is accomplished by setting 
Tc at the optimum value for the larger values of T, thus maximizing Pj) for the star 
magnitude to be mapped. As r decreases, the optimum Tc increases (see discussion 



(a) Cross-correlation method Is detection 
of l's only. 



N 

(b) Cross -correlation method 2: detection 

of 1 1 s and 0 ' s . 

Figure 7*- Probability of detection as a function of 
number of slits and false-alarm probability. 

A = 100; r = 0 . 7 . 
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(a) Cross-correlation method 1: (b) Cross-correlation method 2: 

detection of l's only. detection of I's and O's. 

Figure 8.- Probability of detection as a function of signal-to-noise ratio and 
number of slits. A = 100; = -^”5. 


of fig. 4), and thus Pj) decreases not only because of decreasing T, but also because 
T C is not at its optimum value at the lower r. 

In all cases considered, maximum Pj) occurred for one slit. Thus, use of 
multiple slits did not improve Pp> as was anticipated at the beginning of this study. 
Increasing the number of slits decreases the overall signal-to-noise ratio r/N involved 
in detecting 1 bits, so that should be reduced if Tg is chosen to maintain some 

desired value of Pol- This same effect is carried over to correlation detection, and the 
advantage of increasing code length at fixed Pp\A.,G is a PP aren tly more than offset by 
the degrading effect on detection of l's. 


CONSIDERATIONS FOR PRACTICAL SYSTEMS 


The results of the preceding section apply only to an idealized situation which is not 
possible in practice. In this section the effects of code construction, synchronization 
limitations, and photomultiplier tube noise are studied. Equations are derived which 
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permit the inclusion of these parameters in the analysis of the bit and cross-correlation 
detection processes. 


Code Construction 

In the simplified theoretical analysis the probability of detection was calculated 
according to assumption 3. That is, either the entire code (not necessarily correct) from 
the presence of a star was properly placed in the correlator, or only samples of back- 
ground noise were in the correlator. Thus only the code length, not the arrangement, 
affected the resulting statistics. 

The coding arrangement does affect the probability of false alarm as the detected 
code of a star is shifted into and out of the correlation process. A false alarm during 
this process will reflect directly as an error or uncertainty in the indicated star position. 
Consequently, it is desirable to choose an arrangement that minimizes the rms position 
error. This requires that the probability of a detection be calculated for each code off- 
set in the correlator. The code offset is defined as the number of positions in the cor- 
relator before (-) or after (+) the situation where the received code is properly placed for 
a direct comparison with the known code, that is, the correct star position. Calculating 
the probability of detection for each code offset can be a lengthy and complicated process 
since all possible combinations that can result in a false detection must be taken into 
account. The. results of such a calculation from equation (9) for two codes are presented 
in table I. 

An arbitrary code length of six bits was chosen in keeping with assumption 7. The 
criterion for selection of the two codes 110010 and 110100 was the minimum probability of 
false synchronization for a six-bit code (ref. 6). The codes are compared for 


TABLE I.- ROOT-MEAN-SQUARE POSITION-ERROR CALCULATIONS FOR CODE COMPARISON 


Received 





Code offset in 
correlator, 

6 

Probability of detection for 
various code offset positions, 

P D,6 

Root -mean- square position- error calculations 

digit . . . 

1 

0 

1 

0 

Code 110010 

Code 110100 

Desired 
digit . . . 

1 

0 

0 

1 

Code 

110010 

Code 

110100 

Method 1 

Method 2 

Method 1 

Method 2 

Method 1 

Method 2 

c 

0 

3 

0 

3 

5 

4, 5 

0.9928 X 10-5 

0.9984 x 10-5 

24.82 x 10-5 

24.96 x 10-5 

40.70 x 10-5 

40.93 x 10-5 

— ' c 
w o 

0 

2 

1 

3 

-5, 2 

-5 

.9928 

.1260 

28.79 

3.65 

24.82 

3.15 

%Tl 
XJ rt 

1 

2 

1 

2 

1, -3, -4 

1, -3 

292.76 

36.171 

7611.76 

940.45 

2927.60 

361.71 

sl! 

0 

1 

2 

3 

-2 

-4 

.9928 

.0654 

3.97 

.26 

15.88 

1.05 

0) 0) 
JD > 

1 

1 

2 

2 

-1 

-2, -1 

292.76 

7.2083 

292.76 

7.21 

1463.80 

36.04 

g 2 

1 

3 

0 

2 

3, 4 

2, 3 

292.76 

301.1739 

7319.00 

7529.35 

3805.88 

3915.26 

fc * 

3 

3 

0 

0 

0 

0 

P D,1 

P D,2 

0 

0 

0 

0 








Mean square, = 

15281.10 X 10~ 5 

8505.88 X 10“5 

8278.68 x 10" 5 

4358.14 X 10“ 5 









£2 = 

0.1528110 

0.0850588 

0.0827868 

0.0435814 








Root mean square, 6 = 

0.3909104 

0.291649 

0.287727 

0.2087616 
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Pp^ ~ 10 -5 , X = 100, and T = 0.7. From the computer program results, PgQ, ^11’ 
and T0 are known and are applied to each of the possible combinations of received and 
known bit patterns for detection probability. In the following example the upper code is 
the known code and the lower code represents a received star code preceded and followed 
by 0's which are represented by x's: 


x 


110010 


xxx 


110 


OlOxx 


In this example the code offset position in the correlator is -3, and as the first four col- 
umns of table I indicate for this case, there is one agreement of l's and two agreements 
of 0's; for perfect agreement, another 1 would be detected as a 0 and two more 0's would 
be detected as l's. A detection would occur whenever any combination of these situations 
would cause Tq to be equaled or exceeded. Once Ppg was calculated for each off- 
set, the rms position error relative to zero offset was calculated from the relation 

2N-1 

* 2 = 2 62p D,e < 21 > 

6=1-2N 


where 

6 offset of detected code in correlator 

Pd g probability of detection for offset 6 

e rms position error relative to correct star position 

Table I shows that the lowest e occurs with the 110100 code and method 2 cross- 
correlation detection. No matter which correlation method was used, the same code was 
superior. Also, for either code, method 2 was superior. The reason for this can be 
seen in figure 9, where Pj^g is shown as a function of 6 for each of the situations. 

The dashed line indicates the level of Pp^ w h en n0 star code elements are in the cor- 
relator. Note that for method 2, where both 0's and l's are cross correlated, P D g can 
fall below the dashed line, thus decreasing the rms position error. Because of this char- 
acteristic, errors due to shifting detected star codes into the correlator are less likely 
for method 2 cross-correlation detection. 

Since the calculation of P D.6 is rather lengthy, it is convenient to make a quali- 
tative comparison of codes by plotting the curves shown in figure 10. Here it is assumed 
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Code 110010 


Code 110100 




Code offset in correlator Code offset in correlator 

(a) Cross-correlation method 1. 

Code 110010 Code 110100 




Code offset in correlator Code offset in correlator 

(b) Cross-correlation method 2. 

Figure 9- - Probability-of-detection curves for example codes. The dashed 
line indicates the level of Pp^ when no star code elements are in 
the correlator. 


that the received code is perfectly detected and is preceded and followed by all 0's. For 
method 1 the total 1 bit agreements are plotted against code position, while for method 2 
the total 1 plus 0 bit agreements are used. From a comparison of figures 9 and 10, it is 
evident that the corresponding curves have the same shapes. Hence, qualitative code 
comparisons could be made by comparing agreement levels and curve spreads such as 
those plotted in figure 10. 

For the two codes studied, the probability of false detection is less with negative 
code offsets than with positive code offsets for method 2. For those cases, if two or more 
detections occur during one entire code correlation process, the first detection should be 
chosen as the star position. It is important to note that although method 1 cross correla- 
tion detects only l's, the arrangement of l's and 0's still affects the probability-of- 
detection curves, as shown in figure 9. 
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Code 110010 


Code 110100 




Code offset in correlator 


(a) Cross-correlation method 1. 


Code 110010 


Code 110100 




Code offset in correlator 


(b) Cross-correlation method 2. 

Figure 10.- Agreement -level curves for example codes. 


Synchronization Limitations 

Although in the theoretical analysis it is assumed that synchronization is known, in 
the actual system it is impossible to know when to expect a star pulse. However, it is 
possible to know the repetition rate within limited accuracy. This can be calculated from 
the vehicle spin rate u>, which is obtained by any suitable means. The photomultiplier 
output can be sampled several times faster than the calculated bit rate so that a sufficient 
number of integrated samples are obtained within the star pulse width time tg. It is 
desirable to know how much faster than the repetition rate the sampling rate must be, 
within what accuracy the repetition rate must be known, and what effect these factors have 
on the original theoretical prediction. 

Let the following describe a typical code pattern: 
f r repetition bit rate of actual code to be sampled 
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ideal repetition rate of sampling process 


sfj 1 

t'sfr actual repetition rate of sampling process where v is the error factor 

2N total number of bits in the code 


Let the original synchronized code to be sampled be represented by a square wave 
as shown in figure 11(b). For s = 5, figure 11(c) shows the results of a perfectly sam- 
pled code, that is, both f r and synchronization are known. Figure 11(d) shows three 
possible unsynchronized cases, with uncertain samples indicated by hatching. If synchro- 
nization is not known but f r is, then it is possible that as many as one sample out of 
every five will not have been taken accurately. This condition is independent of code 
length since f r is known, and in general (s - 1) samples per original code bit can be 
utilized. 

Allowing sf r to lie in error will make the sampling process occur too fast or too 
slow according to the error factor v. If every sth (every 5th in fig. 11) sample is com- 
pared with the bits of the known code, then there is a maximum and minimum v for 
which a proper correlation will occur for at least one of the sets of every sth sample. 
Figure 11(e) shows this situation where the first good sample of each bit is the one to be 
correlated. As many as three samples at the end of the sampled code could be in error 
and still permit a proper correlation. Such a situation would occur when < f r . Sim- 
ilarly, figure 11(f) utilizes the last good sample of each bit. Again as many as three 
samples at the end of the sampled code could be in error. This situation would occur 
when i4 r > f r . Since the sampled bits are shifted through the correlator in a serial man- 
ner, bit-by-bit fashion, both correlation situations described above will occur during each 
five correlation processes. Therefore, the proper correlation is guaranteed if the error 

is kept within the allowed limits, which are ± s - seconds for a code 2N bits long. 

s \f r / 

The desired error limitation can now be calculated. From the previous discussion 
the limits of v can be expressed as 

1 5 j, s f < 1 

t s + At t s - At 


1 



sf r 2Ns 



sf r 2NsJ 


1 + 


s - 2 
2Ns 


^ u ^ 


1 


1 - 


s - 2 
2Ns 
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I Mill 


„ . = L_ 

min c < 

i + .^ 

2Ns 

where 

s integral multiple of exact code repetition rate 

t s sample length in time 

At uncertainty in sample time 

If of relations (22) is set equal to 1, then s = 2. Thus a sampling rate of' at 

least twice the code repetition rate is required even if f r is perfectly known. In fig- 
ure 12 is shown as a function of sampling rate s for various code lengths. It is 

readily seen that for a given s, v -m\ n approaches 1 as N increases (the error per- 
mitted in the sampling rate measurement decreases). Thus for a given accuracy limita- 
tion a maximum usable code is established. 



Sampling rate per bit, s 

Figure 12.- Error factor as a function of sampling rate per "bit for various 

code lengths 2N. 

Unfortunately, the synchronization limitations have not only affected the maximum 
allowable code length, but also degraded the probability of detection. The integration 

t k t 

time has been divided by s. In other words, now r = — and A. = - n B (from eq. (2)). 

s s 

It should be recalled that the conclusion from the theoretical analysis was that decreasing 
A decreased the probability of detection. 
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Photomultiplier Tube Noise 


Dark current .- Spontaneous emission of electrons (those that occur when the photo- 
detector is completely shielded from external radiation) within the photomultiplier con- 
stitutes the source of dark current (ref. 7). Since the detection statistics used the param- 
eter of photoelectrons per sample (eq. (2)), the effect of dark current should be analyzed 
by using an equivalent dark- current photoelectron rate. This is 

D = -2# (23) 

eG 

where 

D equivalent average dark-current rate, photoelectrons per sample 

I D average anode dark current, amperes 

e electron charge, coulombs 

G average multiplier gain 

Dark current is not a function of the number of slits N, but is a function of the tube 
design, operating voltages, and environment. Thus dark current is included in the theo- 
retical analysis by adding D to each noise term of equation (2). Obviously it has the 
same effect as background light (although it is not a function of N) and will decrease the 
probability of detection. 

Secondary emission noise . - The use of a discrete photoelectron counter as a detector 
was implicit in the assumption that the electron multiplier gain was the same for each 
photoelectron. Thus, there was no difference between counting photoelectrons at the 
cathode and counting bundles of electrons at the anode. In reality such an assumption is 
not valid. 

Much of the early work with electron multipliers involved investigations of the vari- 
ation of gain due to the statistics of the secondary emission process. These studies 
resulted in the conclusion that the noise was a function of several parameters, including 
the emitter surface and the total number of dynodes (ref. 8). It was also shown that the 
use of the Poisson distribution to describe the gain of each dynode did not account for all 
the additional noise (ref. 9). 

The total variance may be written as 

cr 2 = kT _Z_ b _ kT (j + B ) (24) 

1 y - 1 
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where 


kr variance of conventional shot noise (eq. (1)) 

y average gain per dynode 

b additional noise factor 

1 + B = = ^— b 

y - 1 


The term — was theoretically derived by Morton on the assumption that Poisson 

7-1 

statistics applied to secondary emission (ref. 10). Since this is not adequate, as 
described above, the term b was included to make the equality valid. Winans and 
Pierce used the term 1 + B and from laboratory data believed that it usually lies 
between 1.5 and 3.0. They noted that 1 + B can be determined only by experiment 
(ref. 8). 

Obviously the additional variance will degrade the system performance. For a 
given false-alarm probability, the result will be a significant reduction of the probability 
of detection. 


The effect of limited synchronization knowledge and of photomultiplier noise can be 
included in the statistical analysis through use of the following equations: 

Noise variance, 


(^ + d)(1 + B) 


(25a) 


Signal variance, 


l^ 1 + B ) 

Signal and noise variance, 


(25b) 


§(N + T) + Dj 


(1 +B) 


(25c) 


CONCLUDING REMARKS 


An investigation was made to determine the optimum detection scheme for a star- 
field mapping system that uses coded detection resulting from starlight shining through 
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specially arranged multiple slits of a reticle. It was determined that an optimum scheme 
should be based on the Neyman-Pearson criterion (fixed false-alarm probability) rather 
than the Bayes criterion. With threshold levels established by the Neyman-Pearson cri- 
terion, computer solutions of equations derived from the simplified theory showed that 
the probability of detecting a star could not be improved by using multiple slits in the 
modulating reticle. However, interpretation of the data showed that the use of multiple 
slits improved the system's ability to reject undesirable lower intensity stars. The num- 
ber of slits for optimum rejection is not explicit, since it involves a trade off between the 
probability of detection desired for the stars to be rejected and the minimum allowable 
probability of detection for the lower intensity stars to be mapped. 

A study of code construction showed that some codes are superior to others in 
terms of rms star-position error. It was demonstrated that method 2 (both 0 and 1 cor- 
relation) offered better position error characteristics than method 1 (only 1 correlation) 
for the cases analyzed. Since analysis showed that method 2 was also superior in rejec- 
tion properties, it can be said that method 2 cross-correlation is superior to method 1 
cross -correlation. Also, comparative results indicate that codes can be studied quali- 
tatively with a simple technique of plotting the 1 and 0 agreements for various code 
offsets. 

Analysis of the synchronization problem showed that it was possible to sample at a 
multiple of the expected repetition rate even if the rate was not known perfectly. The 
repetition rate error was shown to determine the necessary sampling rate and the maxi- 
mum code length that could be utilized. The effect of multiple sampling was to decrease 
the resulting probability of detection. 

Equations that were used for the simplified theoretical analysis were modified to 
take into account photomultiplier-tube dark-current and secondary-emission noise. The 
use of these equations in a new analysis, with experimentally determined values for the 
noise parameters, would result in lower probabilities of detection. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 25, 1970. 
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APPENDIX 


JUSTIFICATION FOR COUNTING PHOTOELECTRONS 


The purpose of this appendix is to justify the simplification of the detection process 
to that of counting photoelectrons. The following definitions apply for this appendix only: 

A effective area of optics 

c speed of light 

E(A) energy, joules per photon 

H(A) transmission of optics as a function of A 

h Planck's constant 

P total average power 

P^ average power spectral density 

p average number of photoelectrons per second 

R^ spectral radiance of a star 

r(A) average number of photons per second 

6(A) average quantum efficiency of photocathode as a function of A 

A wavelength 

The energy associated with one photon is 

E(A)=^ (Al) 


Thus the average arrival rate of photons is related to the average power by 


r(A) = 


P 

E(A) 



(A2) 
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APPENDIX — Concluded 


Since the spectral radiance of a star is expressed as an average power spectral density, 
the average number of photoelectrons per second for a given dX becomes 


dp = 


(PX <**) 
he 


(A3) 


where 


P X = R X AH(X)6(X) (A4) 

The total average number of photoelectrons emitted from the photocathode is 


P = dp = A j* R x H(X)6(X)dX 


(A 5) 


where p — k in the theoretical analysis (eq. (1)). 

Note that although <5(X) represents an average quantum efficiency at a given X, 
the use of Poisson statistics is not altered. This was shown by Steinberg and LaTourette 
(ref. 11). 
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